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Gamma-ray induced air showers are notable for their lack of muons, compared to hadronic showers. 
Hence, air shower arrays with large underground muon detectors can select a sample greatly enriched 
in photon showers by rejecting showers containing muons. IceCube is sensitive to muons with 
energies above ~500 GeV at the surface, which provides an efficient veto system for hadronic air 
showers with energies above 1 PeV. One year of data from the 40-string IceCube configuration was 
used to perform a search for point sources and a Galactic diffuse signal. No sources were found, 
resulting in a 90% C.L. upper limit on the ratio of gamma rays to cosmic rays of 1.2 x 10 -3 for the 
flux coming from the Galactic Plane region ( —80° < I < —30°; —10° < b < 5°) in the energy range 
1.2 - 6.0 PeV. In the same energy range, point source fluxes with E~ 2 spectra have been excluded at 
a level of (£/TeV) 2 d$/dS ~ 10" 12 - 10" 11 cm" 2 s" 1 TeV" 1 depending on source declination. The 
complete IceCube detector will have a better sensitivity, due to the larger detector size, improved 
reconstruction and vetoing techniques. Preliminary data from the nearly-final IceCube detector 
configuration has been used to estimate the 5 year sensitivity of the full detector. It is found to 
be more than an order of magnitude better, allowing the search for PeV extensions of known TeV 
gamma-ray emitters. 



I. INTRODUCTION 

Gamma-rays are an important tool for studying the 
cosmos; unlike cosmic rays (CRs), they point back to 
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their sources and can identify remote acceleration re- 
gions. Air Cherenkov telescopes have identified numer- 
ous sources of high-energy (E > 1 TeV) gamma-rays (see 
e.g. pQ): within our galaxy, gamma-rays have been ob- 
served coming from supernova remnants (SNRs), pulsar 
wind nebulae (PWNe), binary systems, and the Galactic 
Center. Extra-galactic sources include starburst galaxies 
and Active Galactic Nuclei (AGNs). Surface air-shower 
arrays like Milagro have performed all-sky searches for 
TeV gamma-rays. Although these detectors are less sen- 
sitive to point sources than Air Cherenkov telescopes, 
they have identified several Galactic pointlike and ex- 
tended sources [3J. Interactions of CRs with interstellar 
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matter and radiation in the Galaxy produce a diffuse flux. 
Hadrons interacting with matter produce neutral pions, 
which decay into gamma rays, while CR electrons pro- 
duce gamma-rays via inverse Compton scattering on the 
radiation field. Milagro has measured this diffuse Galac- 
tic flux in the TeV energy range with a median energy 
of 15 TeV and reported an excess in the Cygnus region, 
which might originate from CRs from local sources inter- 
acting with interstellar dust clouds [S]- IceCube's prede- 
cessor AMANDA-II has also looked for TeV photons from 
a giant flare from SGR 1806-20, using fOO GeV muons. 
AMANDA's large muon collection area compensated for 
the small cross-section for photons to produce muons [I] . 

At higher energies, extra-galactic sources are unlikely 
to be visible, because more energetic photons are pre- 
dicted to interact with cosmic microwave background ra- 
diation (CMBR), and with infrared starlight from early 
galaxies, producing e + e~ pairs [5]. At 1 PeV, for exam- 
ple, photon propagation is limited to a range of about 10 
kpc. It is unknown whether Galactic accelerators exist 
that can produce gamma rays of such high energy, but 
an expected flux results from interaction of (extragalac- 
tic) CRs with the interstellar medium (ISM) and dense 
molecular clouds. 

To date, the best statistics on photons with energies 
in the range from ~300 TeV to several PeVs come from 
the Chicago Air Shower Array - Michigan Muon Array 
(CASA-MIA), built at the Dugway Proving Ground in 
Utah. CASA consisted of 1089 scintillation detectors 
placed on a square array with 15 m spacing. MIA con- 
sisted of 1024 scintillation counters buried under about 
3 m of earth, covering an area of 2500 m 2 . It served as a 
muon veto, with a threshold of about 0.8 GeV. 

CASA-MIA set a limit on the fraction of photons in 
the cosmic-ray flux of 10 -4 at energies above 600 TeV 
[5J. The experiment also sets a limit of 2.4 x 10~ 5 on the 
fraction of photons in the CR flux coming from within 5° 
of the galactic disk [7] at 310 TeV. This is near the the- 
oretical expectation due to cosmic-ray interactions with 
the interstellar medium. For a Northern hemisphere site 
like CASA-MIA, Ref. [8] predicts a gamma-ray fraction 
of 2 x 10~ 5 for the average gas column density. 

In this work, we present a new approach for detect- 
ing astrophysical PeV gamma rays, based on data of 
the surface component, IceTop, and the in-ice array of 
IceCube. IceTop measures the electromagnetic compo- 
nent of air showers, while the in-ice array is sensitive 
to muons that penetrate the ice with energies above 
500 GeV. While most CR showers above 1 PeV contain 
many muons above this threshold, only a small fraction of 
PeV gamma-ray showers carry muons that are energetic 
enough to reach the in-ice array. Therefore, gamma-ray 
candidates are selected among muon-poor air showers de- 
tected with IceTop and whose axis is reconstructed as 
passing through the in-ice array. 

This approach of selecting muon-poor showers as 
gamma-ray candidates is fundamentally different from 
the earlier AMANDA-II gamma-ray search described 




FIG. 1. IceCube consists of a ~km 2 surface air shower 
array and 86 strings holding 60 optical modules each, filling 
a physical volume of a km 3 . The region in the center of the 
buried detector is more densely instrumented. See text for 
details. 



above, which was only sensitive to gamma-ray showers 
that do contain high energy muons (> 100 GeV). 

We present a limit on the gamma-ray flux coming from 
the Galactic Plane, based on one year of data with half 
of the IceCube strings and surface stations installed. We 
also discuss the sensitivity of the completed detector. 



II. THE ICECUBE DETECTOR 

IceCube (see Fig. [lj is a particle detector located at 
the geographic South Pole. The in-ice portion consists 
of 86 strings that reach 2450 m deep into the ice. Most 
of the strings are arranged in a hexagonal grid, sepa- 
rated by ~125 m. Each of these strings holds 60 digital 
optical modules (DOMs) separated by ~ 17 m covering 
the range from 1,450 m to 2,450 m depth. Eight strings 
form a denser instrumented area called DeepCore. The 
DOMs detect Cherenkov light produced by downward- 
going muons in cosmic-ray air showers and from charged 
particles produced in neutrino interactions. The data 
used in this analysis was collected in 2008/9, when the 
40 strings shown in Fig. [2] were operational. 

Each DOM is a complete detector system, comprising 
a 25 cm diameter Hamamatsu R7081-02 phototube [9], 
shaping and digitizing electronics |10j . calibration hard- 
ware, plus control electronics and power supply. Most of 
the buried PMTs are run at a gain of 10 7 . Digitization is 
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initiated by a discriminator, with a threshold set to 0.25 
times the typical peak amplitude of a single photoelec- 
tron waveform. Each DOM contains two separate digi- 
tizing systems; the Analog Transient Waveform Digitizer 
(AT WD) records 400 ns of data at 300 Megasamples/s, 
with a 14 bit dynamic range (divided among 3 paral- 
lel channels), while the fast Analog-to-Digital Converter 
(fADC) records 6.4 fis of data at 40 Megasamples/s, with 
10 bits of dynamic range. A system transmits timing 
signals between the surface and each DOM, providing 
timing calibrations across the entire array of about 2 ns 

mm- 

The IceTop surface array [12] is located on the surface 
directly above the in-ice detectors. It consists of 81 sta- 
tions, each consisting of two ice-filled tanks, about 5 m 
apart. For the 2008 data used here, 40 stations were op- 
erational (IC40, see Fig. [2]). Each tank is 1.8 m in diam- 
eter, filled with ice to a depth of about 90 cm. The tanks 
are initially filled with water, and the freezing of the wa- 
ter is controlled to minimize air bubbles and preserve 
the optical clarity of the ice. Each tank is instrumented 
with two DOMs, a high-gain DOM run at a PMT gain 
of 5 x 10 6 , and a low-gain DOM, with a PMT gain of 
5 x 10 5 . The two different gains were chosen to maximize 
dynamic range; the system is quite linear over a range 
from 1 to 10 5 photoelectrons. A station is considered hit 
when a low-gain DOM in one tank fires in coincidence 
with a high-gain DOM in the other; the thresholds are 
set to about 20 photoelectrons. 

When an in-ice DOM is triggered it sends a Local Coin- 
cidence (LC) message to its nearest two neighbors above 
and nearest two neighbors below. If the DOM also re- 
ceives an LC message from one of its neighbors within 1 
fj,s it is in Hard Local Coincidence (HLC). In that case 
the full waveform information of both the ATWD and 
fADC chip is stored. For IC40 and earlier configurations, 
isolated or Soft Local Coincidence (SLC) hits were dis- 
carded. In newer configurations, the SLC hits are stored 
albeit with limited information. Keeping the full wave- 
forms would require too much bandwidth, since the rate 
of isolated hits per DOM due to noise is ^500 Hz. In- 
stead, only the three fADC bins with the highest values 
and their hit times are stored. 

In Sec. |IV| we present a gamma- ray analysis of the IC40 
dataset. In Sec. [VI] we study the expected sensitivity of 
the full IceCube detector, and discuss how the inclusion 
of SLC hits increases the background rejection. 



III. DETECTION PRINCIPLE 

We create a sample of gamma-ray candidates by select- 
ing air showers that have been successfully reconstructed 
by IceTop and have a shower axis that passes through 
the IceCube instrumented volume. The geometry limits 
this sample to showers having a maximum zenith angle of 
^30 degrees. Since IceCube is located at the geographic 
South Pole, the Field of View (FOV) is roughly con- 
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FIG. 2. Map of location of all 86 strings of the completed 
IceCube detector. The blue dots represent the 40 string con- 
figuration that is used for this analysis. At surface level each 
of these 40 strings is complemented by an IceTop station con- 
sisting of two tanks. The large (red) circles indicate the 'inner 
strings' of the IC40 configuration. 



strained to the declination range of —60 to —90 degrees, 
as shown in Fig. [3] This FOV includes the Magellanic 
Clouds and part of the Galactic Plane. Gamma-rays at 
PeV energies are strongly attenuated over extra-galactic 
distances, thus limiting the observable sources to those 
localized in our galaxy. At distances of ~ 50 kpc and 
~60 kpc, the PeV gamma-ray flux from the Large and 
Small Magellanic Cloud is suppressed by several orders 
of magnitude. 

The contours in the background of Fig. [3] are the in- 
tegrated neutral atomic hydrogen (HI) column densities 
under the assumption of optical transparency based on 
data from the Leiden/ Argentine/Bonn survey [H] . These 
densities are not incorporated into the analysis and are 
only plotted to indicate the Galactic Plane. We do ex- 
pect, however, that gamma-ray sources are correlated 
with the HI column density. Firstly, Galactic CR acceler- 
ators are more abundant in the high density regions of the 
Galaxy. Secondly, the gamma-ray flux of (extra-)galactic 
CRs interacting with the ISM naturally correlates with 
the column density. However, it has to be noticed that 
this correlation is not linear, because of the attenuation of 
gamma-rays over a 10 kpc distance scale. Furthermore, 
the column densities do not include molecular hydrogen 
which can also act as a target for CRs. 

The gamma-ray candidate events are searched for in 
a background of CR showers that have exceptionally few 
muons or are directionally misreconstructed. In the latter 
case the muon bundle reaches kilometers deep into the ice 
but misses the instrumented volume. This background 
is hard to predict with Monte Carlo (MC) simulations. 
Cosmic-ray showers at PeV energies and with a low num- 
ber of energetic muons are rare. For example, at 1 PeV 
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FIG. 3. Contours of integrated neutral atomic hydrogen 
(HI) column densities [14], in Galactic coordinates (flat pro- 
jection). The blue circle indicates the gamma-ray FOV for 
IceCube in the present IC40 analysis. The red rectangle in- 
dicate the regions for which CASA-MIA [B] has set an upper 
limit on the Galactic diffuse photon flux in the 100 TeV — 
1 PeV energy range. IceCube's FOV is smaller but covers 
a different part of the Galactic Plane, close to the Galactic 
Center. 

less than 0.1% of the simulated showers contain no muons 
with an energy above 500 GeV, approximately what is 
needed to reach the detector in the deep ice when trav- 
eling vertically. Determining how many hadronic show- 
ers produce a signal in a buried DOM would require an 
enormous amount of MC data to reach sufficient statis- 
tics, plus very strong control of the systematic uncer- 
tainties due to muon production, propagation of muons 
and Cherenkov photons through the ice, and the abso- 
lute detector efficiencies. It would also have to be able to 
accurately predict the errors in air shower reconstruction 
parameters. For example, this analysis is very sensitive 
to the tails of the distribution of the error on the an- 
gular reconstruction of IceTop. Even MC sets that are 
large enough to populate these tails are not expected to 
properly describe them. 

To avoid these issues, we determine the background di- 
rectly from data. As a result, we are not able to measure 
a possible isotropic contribution to the gamma-ray flux, 
because these gamma-rays would be regarded as back- 
ground. Instead, we search for localized excesses in the 
gamma-ray flux. 

We search for a correlation of the arrival directions of 
the candidate events with the Galactic Plane, and scan 
for point sources. The acceptance of IceTop-40 is a com- 
plex function of azimuth and zenith due to its elongated 
shape and the requirement that the axis of the detected 
shower passes through IC40 (with the same elongated 
shape). However, since the arrival time is random (there 
are no systematic gaps in detector uptime w.r.t. side- 
real time) the reconstructed right ascension (RA) of an 
isotropic flux of CR showers is uniform. The correct dec- 
lination distribution of the background is very sensitive to 



the number of background showers introduced by errors 
in the IceTop angular reconstruction of the air shower 
as a function of the zenith angle, and is taken to be un- 
known. However, the flat distribution of background over 
RA is enough to allow for a search for gamma-ray sources. 

Recently, IceCube found an anisotropy in the arrival 
direction of CRs on the southern hemisphere [15] . These 
deviations with RA have been established for samples 
of CRs with median energies of 20 TeV and 400 TeV. 
The two energy ranges show a very different shape of the 
anisotropy, but the level of the fractional variations in 
flux is at a part-per-mil level for both |16j . An anisotropy 
with comparable magnitude in the PeV energy range is 
too small to affect this analysis (the IC40 final sample 
contains 268 events). 

IV. IC40 ANALYSIS 

A. Event selection 

Between April 5 2008 and May 20 2009, IceCube took 
data with a configuration of 40 strings and 40 surface 
stations (IC40) , using several trigger conditions based on 
different signal topologies. This analysis uses the 8 sta- 
tion surface trigger, which requires a signal above thresh- 
old in both tanks of at least 8 IceTop stations. An ad- 
ditional signal in IceCube is not required for this trig- 
ger, but all HLC hits in the deep detector within a time 
window of 10 fj,s before and after the surface trigger are 
recorded. 

The air shower parameters are reconstructed from 
the IceTop hits with a series of likelihood maximization 
methods. The core position is found by fitting the lateral 
distribution of the signal, using 



where S is the signal strength, r is the distance to the 
shower axis, S rc { is the fitted signal strength at the ref- 
erence distance i? ro f = 125 m, j3 is the slope parameter 
reflecting the shower age, and k = 0.303 is a constant de- 
termined from simulation 13, 23 . Signal times are used 
to find the arrival direction of the air shower. The time 
delay due to the shape of the shower plane is described by 
the sum of a Gaussian function and a parabola, both cen- 
tered at the shower core, which yields a resolution of 1.5° 
for IC40. The relationship between the reconstructed en- 
ergy E Teco and S Te f is based on MC simulations for proton 
showers and depends on the zenith angle. 

IceCube data is processed in different stages: in the 
first two levels the raw data is calibrated and filtered, 
and various fitting algorithms are applied, of which only 
the IceTop reconstruction described above is used in this 
analysis. In the selection of photon shower candidates 
from the data sample we distinguish two more steps: level 
three (L3) and level four (L4) . Level three includes all the 



conditions on reconstruction quality, geometry and en- 
ergy that make no distinction between gamma-ray show- 
ers and CR showers. The L4 cut is designed to separate 
gamma rays from CRs. 

Two parameters are used to constrain the geometry 
and ensure the shower axis passes through the instru- 
mented volume of IceCube. The IceTop containment pa- 
rameter Cit is a measure for how centralized the core 
location is in IceTop. When the core is exactly in the 
center of the array Cit = 0, while Cit = 1 means that 
it is exactly on the edge of the array. More precisely, 
Cit = x means that the core would have been on the 
edge of the array if the array would be x times its actual 
size. The string distance parameter d str is the distance 
between the point where the shower reaches the depth 
of the first level of DOMs and the closest inner string. 
Inner string, in this sense, means a string which is not 
on the border of the detector configuration. IC40 has 17 
inner strings (see Fig.j2|. The L3 cuts are: 

• Quality cut on lateral distribution fit: 1.55 < (3 < 
4.95 (cf. Eq.[T]) 

• Geometry cut: Cit < C max 

• Geometry cut: d s tr < d max 

• Energy cut: E ICCO > E min 

The energy and geometry cuts are optimized in a later 
stage (Sec. llVC]). 



The L4 stage imposes only one extra criterion to the 
event: there should be no HLC hits in IceCube. This 
removes most of the CR showers, if E m - m is chosen suffi- 
ciently high. The remaining background consists of CR 
showers with low muon content and mis-reconstructed 
showers, the high energy muons from which do not actu- 
ally pass through IceCube. 

The event sample after the L4 cut might be dominated 
by remaining background but it can be used to set an 
upper limit to the number of gamma-rays in the sample. 
Since the event sample after L3 cuts is certainly domi- 
nated by CRs, the ratio between the number of events 
after L4 and L3 cuts can be used to calculate an upper 
limit to the ratio of gamma-ray-to-CR showers. 

The remaining set of candidate gamma ray events 
is tested for a correlation with the Galactic Plane 
(Sec. IV C I and the presence of point like sources 
(Sec. |IVD ). First, the results of simulations are pre- 



sented, which provide several quantities needed for sen- 
sitivity calculations. 



B. Simulation 

Although we determine the background from data only, 
simulations are needed to investigate systematic differ- 
ences in the detector response to gamma-ray showers and 
cosmic-ray showers. More specifically, we are interested 
in the energy reconstruction of gamma-ray showers, the 
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FIG. 4. Ratio between reconstructed and true energy of sim- 
ulated gamma-ray showers as a function of their true energy. 
At low energies the overestimation of the gamma-ray energy 
is largely due to a bias effect of the eight-station filter. At 
higher energies, this overestimation decreases. 



fraction of gamma-ray showers that is rejected by the 
muon veto system, and a possible difference in effective 
detector area between both types of showers. 

Gamma-ray and proton showers are simulated with 
CORSIKA v6.900, using the interaction models FLUKA 
2008 and SYBILL 2.1 for low and high energy hadronic 
interactions, respectively. For both primaries 10,000 
showers are generated within an energy range of 800 TeV 
to 3 PeV with a E^ 1 spectrum. Because the shower axes 
are required to pass through IceCube, the zenith angle 
is restricted to a maximum of 35 degrees. The observa- 
tion altitude for IceTop is 2835 m. Atmospheric model 
MSIS-90-E is used, which is South Pole atmosphere for 
July 01, 1997. Seasonal variations in the event rate are 
of the order of 10% Q2]. 

The detector simulation is done with the IceTray soft- 
ware package. Each simulated shower is fed into the de- 
tector simulation ten times with different core positions 
and azimuthal arrival direction, for a total of 100,000 
events for both gamma rays and protons. This resam- 
pling of showers is a useful technique for increasing the 
statistics when examining quantities like the resolution 
of the energy reconstruction of IceTop. However, it can- 
not be used for quantities with large shower-to-shower 
variations, such as the number of high energy muons. 

The energy of gamma-ray showers is overestimated by 
the reconstruction algorithm. Fig. [4] shows the distri- 
bution of the logarithm of the ratio between the recon- 
structed and true primary energy as function of true en- 
ergy, weighted to a _E~ 2 ' 7 spectrum. There are two rea- 
sons for the energy offset. First, there is a selection ef- 
fect of the eight-station filter, which has a bias (below a 
few PeV) towards showers that produce relatively large 
signals in the IceTop tanks. This also affects the recon- 
structed energy of CR showers. At higher energies, the 
offset decreases but the reconstructed energy of gamma- 
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FIG. 5. Distribution of true energy of gamma-ray (red, 
solid) and proton (blue, dotted) showers for an energy cut at 
E leco > 1.4 PeV, weighted to a E~ 2 ' 7 spectrum. 



ray showers is still slightly overestimated because the en- 
ergy calibration of IceTop is performed with respect to 
proton showers. 

Figure [5] shows the distribution of true energies of 
gamma-ray and proton showers for the energy cut 
E leco > 1.4 PeV (which will be adopted in Sec. IV C I. 
After this cut, 95% of the gamma-ray showers have a true 
energy above 1.2 PeV, while 95% of the proton showers 
have an energy above 1.3 PeV. 

The fraction of gamma-ray showers that is falsely re- 
jected because the showers contain muons that produce 
a signal in IceCube is found by applying the cuts to the 
MC simulation. 



After applying the L3 cuts (defined in Sec. IV C I to the 
simulated gamma-ray sample there are 737 events left in 
the sample, of which 121 produce a signal in IceCube. 
Taking into account an energy weighting of E~ 2 , this 
corresponds to 16%. 

Showers that have no energetic muons can still be re- 
jected if an unrelated signal is detected by IceCube. This 
could be caused by noise hits or unrelated muon tracks 
that fall inside the time window. This noise rate is de- 
termined directly from the data and leads to 14% signal 
loss. The total fraction of gamma-ray showers that is 
falsely rejected is therefore 28%. 

Finally, because the composition, shower size and evo- 
lution of gamma-ray and CR showers are different, one 
might expect a difference in the number of triggered sta- 
tions and the quality of the reconstruction, which could 
lead to different effective areas. Such an effect would be 
of importance when calculating the relative contribution 
of gamma-rays to the total received flux. We compare 
the effective area for gamma rays and protons by count- 
ing the number of events that are present at L3. We 
compensate for the energy reconstruction offset by re- 
ducing the reconstructed gamma-ray energy by a factor 
1.16. The ratio of the effective area for gamma-rays to 
that for protons is then found to be 0.99. 



It should be emphasized that we do not use the sim- 
ulation to determine the number of muons and their en- 
ergy distribution from CR showers. This would require a 
simulation set that includes heavier nuclei instead of pro- 
tons only. Moreover, various hadronic interaction models 
generate significantly different muon fluxes [18]. Instead, 
this analysis estimates the rate of CR showers that do 
not trigger IceCube using the data itself. 



C. Galactic Plane test 

The IC40 data set consists of 368 days of combined 
IceCube and IceTop measurements. The data from Au- 
gust 2009 is used as a burn sample, which means that it 
is used to tune the parameters of the analysis. After this 
tuning the burn sample is discarded and the remaining 
data is used for the analysis. 

IceCube is sensitive to gamma rays above 1 PeV. Ear- 
lier searches by CASA-MIA in a slightly lower energy 
range (100 TeV - 1 PeV) with better sensitivity have not 
established a correlation of gamma-rays with the Galac- 
tic Plane (see Fig. [l4| . For a Galactic diffuse flux below 
the CASA-MIA limit [B] no gamma rays are expected in 
the IC40 burn sample. However, IceCube observes a dif- 
ferent part of the Galactic Plane (see Fig. [3]), close to the 
Galactic Center, so the possibility exists that previously 
undetected sources or local enhancements in CR and dust 
densities create an increase in the flux from that part of 
the sky. 

In order to find a possible correlation of candidate 
events in our burn sample to the Galactic Plane, different 
sets of L3 cut parameters are applied to find a set that 
produces the most significant correlation. Afterwards, 
the cut parameters are fixed and the burn sample is dis- 
carded. The fixed cut parameters are then used in the 
analysis of the rest of the IC40 data set to test whether 
the correlation is still present. Note that these cuts are 
applied at L3, so they affect both the event samples af- 
ter L3 and L4 cuts. This is important because the ratio 
between the number of events after L4 and L3 cuts is 
used to calculate the ratio between gamma ray and CR 
showers. 

There are three cut parameters that are optimized by 
using the burn sample: i? m i n , C max and d max . This is 
done by scanning through all combination of parameters 
within the following range: 600 TeV < E min < 2 PeV 
with steps of 100 TeV, 0.5 < C max < 1.0 with steps of 
0.1, and 50 m < d max < 90 m with steps of 10 m. For 
each combination, the number of events Ns in the burn 
sample after the L4 cut, located in the source region, is 
counted. The source region is defined as within 10 de- 
grees of the Galactic Plane. Then, the data set is scram- 
bled multiple times by randomizing the RA of each data 
point. For each scrambled data set the number of events 
in the source region is again counted. The best combi- 
nation of cut parameters is the set which has the lowest 
fraction of scrambled data sets for which this number is 



equal to or higher than Ns- 

The result of the scan is given in the three panels of 
Fig. [6j For each cut parameter the fraction of scrambled 
data sets that has a number of events in the source region 
equal to or exceeding the amount in the original data 
set is plotted for different cut values. For each plot the 
values of the other two cut parameters are kept constant 
at their optimal value. The actual search is done in three 
dimensions. The ratio is lowest for E > 1.4 PeV, Cn < 
0.8 and d str < 60 m. With this combination of cuts only 
0.011% of the scrambled data sets produce an equal or 
higher number of events in the source region. 

Note that while this procedure of optimizing cuts 
should be effective in the presence of sufficient signal, 
the small fraction obtained here and its erratic behavior 
with changing cut values are consistent with fluctuations 
of the background CR distribution alone, given the large 
number of possible cut combinations which were scanned. 
Nonetheless, the cut parameter values found with this 
procedure seem very reasonable (similar values are found 
with an alternative method, see Sec. VI). 



The optimized cuts are applied to the complete IC40 
data set minus the burn sample. There are 268 candi- 
date events of which 28 are located in the source region. 
Figure [7] is a map of the sky showing all 268 events. The 
colors in the background indicate the integrated HI col- 
umn densities, cf. Fig. [3] (see discussion Sec. III). These 
are meant to guide the eye and are not part of the anal- 
ysis. 

The significance of the correlation with the Galactic 
Plane is tested by producing data sets with scrambled 
RA. An equal or higher number of source region events 
is found in 21% of the scrambled data sets, corresponding 
to a non-significant excess of +0.9<7. 

We follow the procedure of Feldman & Cousins [20] to 
construct an upper limit for the ratio of gamma rays to 
CRs. The background is determined by selecting a range 
of RA that does not contain the source region. Within 
this range the data points are scrambled multiple times 
and for each scrambled set the number of events in a pre- 
defined region of the same shape and size as the source 
region is counted. This yields a mean background of 
24.13 events for the source region. Using a 90% confi- 
dence interval, the upper limit on the number of excess 
gamma rays from this region is 14. 

Since 28% of gamma-ray showers are rejected by the 
veto from the buried detector, the maximum number of 
excess gamma-rays from the Galactic Plane is 14/0.72 = 
19.4. From Fig. [5j it is known that the energy cut cor- 
responds to a threshold of 1.2 PeV for gamma-rays, and 
1.3 PeV for protons. Given that at L3 the sample is dom- 
inated by CR showers, and assuming a CR and gamma- 
ray power law of 7 = —2.7, a 90% C.L. upper limit of 
1.2 x 10~ 3 on the ratio of gamma-ray showers to CR 
showers in the source region can be derived in the en- 
ergy range 1.2 - 6.0 PeV. The upper bound of 6.0 PeV is 
the value for which 90% of the events are inside the en- 
ergy range. This value falls outside the range for which 



gamma-ray showers were simulated. However, there is 
no indication that the energy relation plotted in Fig. [4] 
behaves erratically above 3 PeV. 

This is a limit on the average excess of the ratio of 
gamma-rays to CRs in the source region with respect to 
the rest of the sky, i.e. a limit on the Galactic component 
of the total gamma-ray flux. A possible isotropic compo- 
nent is not included. Systematic uncertainties lead to a 
18% variation of the upper limit, as determined in Sec.|v| 



D. Unbinned point source search 

An additional search for point-like sources tests the 
possibility that a single source dominates the PeV 
gamma-ray sky. This source does not necessarily lie close 
to the Galactic Plane. An unbinned point source search 
is performed on the sky within the declination range of 
-85° to -60°, using a method that follows [21. The 
region within 5 degrees around the zenith is omitted, be- 
cause the method relies on scrambled data sets that are 
produced by randomizing the RA of the events. Close 
to the zenith this randomization scheme fails due to the 
small number of events. 

The data is described by an unknown amount of signal 
events on top of a flat distribution of background events. 
In an unbinned search, a dense grid of points in the sky 
is scanned. For each point a maximum likelihood fit is 
performed for the relative contribution of source events 
over background events. 

For a particular event i the probability density function 
(PDF) is given by 



Pi(ns) = -^rSi 



1 



N J 1 



(2) 



where ns is the number of events that is associated to 
the source, Bi is the background PDF, and 



Si 



1 



2ira 2 



exp 



(3) 



is the two-dimensional Gaussian source PDF, in which 
A "J is the space angle between the event and the source 
test location, and a = 1.5° is the angular resolution of 
IceTop. The background PDF Bi is only dependent on 
the zenith angle, and is derived from the zenith distri- 
bution of the data. For each point in the sky there is a 
likelihood function 



L{n s ) = UiPiins), 
and associated test statistic 

A = -2(log(L(0))-log(L(n s )) ! 



(4) 



(5) 



which is maximized for ns- In the optimization proce- 
dure, ns is allowed to have a negative value, which math- 
ematically corresponds to a local flux deficit. 

The procedure is similar to the search method for the 
neutrino point sources with IceCube |22j , except that the 
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FIG. 6. Optimization scans for cut parameters -E m i n , d max , and C ma x- The fraction of scrambled data sets that perform equal 
or better than the real data set is plotted against cut value. For each plot, the other 2 parameters are kept constant at their 
optimal value. The actual scan was done three dimensionally. For each plot, the shaded region indicates the parameter space 
that is excluded by the optimized cuts. 
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FIG. 7. Equatorial map of the 268 candidate gamma-ray events of the IC40 data set superimposed on HI column densities 
based on |14| . The dotted black curve encloses the source region, denned as within 10 degrees of the Galactic Plane. 



source and background PDF do not contain an energy 
term. Because the range of energies in the event sam- 
ple is relatively small (90% of the events have an energy 
between 1 and 6 PeV), an energy PDF is unlikely to im- 
prove the sensitivity. 

Figure [8] displays a map of the sky with declination 
between —85° and —60° showing the events in this region 
and contours of the test statistic A. The maximum value 
is A = 2.1 at 5 = -65.4° and RA= 28.7°, corresponding 
to a fit of n s = 3.5 signal events. The overall significance 
of this value for A is found by producing 10, 000 scrambled 
data sets by randomizing the RA of each event. Figure [9] 
shows the distribution of A associated to the hottest spot 
in each scrambled data set. The median test statistic 
value for the hottest spot in the scrambled data sets is 
A = 2.7, so the actual data set is consistent with a flat 
background. 

Upper limits on the gamma-ray flux can be derived 
for each point in the sky by assuming that all events are 
gamma-rays. Since many events are in fact muon-poor 
or misreconstructed showers, this leads to a conservative 
upper limit. Because the acceptance of IC40 decreases as 
a function of declination (see Fig. 13), the limit is more 



constraining at lower declination. Figure 10 is a sky map 



of the 90% C.L. upper limit in the energy range E — 
1.2 — 6.0 PeV for E~ 2 source spectra. Point source fluxes 
are excluded at a level of (£7TeV) 2 d$/d£: - 10~ 12 - 
10~ n cm _2 s _1 TcV _1 depending on source declination. 
Corrections for signal efficiency and detector noise are 
taken into account. Systematic uncertainties lead to a 
18% variation of the upper limit. 



V. SYSTEMATIC ERRORS 

Since this analysis derives the background from data, 
the systematic uncertainties due to the background esti- 
mation are small. The previously discussed cosmic-ray 
anisotropy measurement (see Sec. Ill and |16j ) is too 
small to have an impact on this analysis. Since there 
are no systematic gaps in detector uptime with respect 
to sidereal time-of-day in our sample, the coverage of RA 
is homogeneous. 

Therefore, we focus on the systematic uncertainties in 
the signal efficiency, due to uncertainties in the surface 
detector sensitivity, and in the muon production rate for 
photon showers. 

The uncertainty in the surface detector sensitivity is 
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FIG. 9. Distribution of the largest value of A observed in 
each scrambled data set. The red dotted line indicated the 
value for A that corresponds to the hottest spot in the data. 



two different photoproduction models, and finds (for wa- 
ter with a similar atomic number and mass number as 
air), a difference that rises from about 20% at 10 TeV to 
60% at 1 PeV. The bulk of the particles in the shower are 
at lower energies, so we adopt a 20% uncertainty on the 
muon production rate via photoproduction. In addition, 
there is also a contribution of muon pair creation. To 
reach the in-ice DOMs, muons need at least 500 GeV. 
At 1 TeV, the fractional contribution of muon pair cre- 
ation is ~ 10% |25j . Since muon pair production is not 
included in SYBILL 2.1, we arrive at 30% uncertainty in 
total muon production rate. This uncertainty is applied 
to the 16% of photon showers that are lost because they 
contain muons for a final 4.8% uncertainty in sensitivity 
due to the unknown muon production cross-section. 

We add the uncertainties due to detector response and 
muon production in quadrature, and arrive at an overall 
18 % uncertainty in sensitivity. 



studied in Ref. [33]; Table 2 there gives the uncertain- 
ties for hadronic showers as a function of shower energy 
and zenith angle. Although there are differences be- 
tween hadronic and electromagnetic showers, most fac- 
tors that contribute to this figure apply to both types 
of showers. Strongly contributing factors include atmo- 
spheric fluctuations, calibration stability and uncertain- 
ties in response of detector electronics (PMT saturation 
and droop). The contribution from the uncertainty in 
modeling the hadronic interaction is clearly different for 
electromagnetic showers, and is discussed below. For 
E < 10 PeV, and zenith angle less than 30°, there is 
a 6.0% systematic uncertainty in energy, and a 3.5% sys- 
tematic uncertainty in flux. For an E~ 2,7 spectrum, a 
6.0% uncertainty in energy translates into a 17.0% un- 
certainty in flux, or, adding in quadrature, 17.4% flux 
uncertainty. 

The uncertainty in the muon production from hadronic 
showers emerges from theoretical uncertainties. It de- 
pends on the hadronic photoproduction and electropro- 
duction cross-sections for energies between 10 TeV and 
6 PeV. Figure 1 of [H] compares two cross-sections from 



VI. ICECUBE 5- YEAR SENSITIVITY 

The sensitivity of the full IceCube detector to a 
gamma-ray flux from the Galactic Plane benefits from 
multiple improvements that can be made with respect 
to the analysis presented above. In this section we use 
preliminary data from the IC79 configuration (79 strings, 
73 surface stations, 2010/2011) to estimate the sensitiv- 
ity that the full IceCube detector can reach in 5 years. 
Since the full detector (IC86: 86 strings and 81 surface 
stations) is slightly larger than the IC79 configuration, 
the predicted sensitivity will be slightly underestimated. 
Also, the new cuts proposed below are not yet optimized, 
as this would require the actual IC86 data set. 



A. Air shower reconstruction 

This analysis is very sensitive to the quality of the core 
reconstruction. If the shower core is not reconstructed 
accurately, a muon bundle that passes outside the in-ice 
array might be incorrectly assumed to be aimed at the 
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FIG. 10. Sky map of 90% C.L. upper limits on point source flux (_B/TeV) 2 d$/d_E in cm 2 s 1 TeV 1 for E 2 source spectra 



in the energy range E = 1.2 
largest. 



6.0 PeV. The limit is typically more constraining at low declinations where the effective area is 



detector. Because of the absence of a signal in the in-ice 
DOMs, the event is then misinterpreted as a gamma-ray 
candidate. A more accurate core reconstruction algo- 
rithm has been developed for IceTop and will improve 
the CR rejection in post-IC40 analyses. In addition, the 
angular resolution of the larger array is improved, in- 
creasing the sensitivity to point sources. 



B. Isolated hits 

The SLC mode (which is available since IC59, see 
Sec. |n| increases the sensitivity to CR showers with low 
muon content. A muon with just enough energy to reach 
IceCube, might not emit enough Cherenkov light to trig- 
ger multiple neighboring DOMs. By tightening the L4 
cut so that no SLC hits are allowed to be present in the 
data, the efficiency with which CR showers can be re- 
jected increases. At the same time, actual gamma-ray 
showers may be rejected in case of a noise hit in a single 
DOM. To keep this chance low, SLC hits only count as 
veto hits if they can be associated to the shower muon 
bundle both spatially and temporally. 

Figure [TT] shows the distribution of isolated hits in the 
complete detector as a function of time relative to the 
arrival time of the air shower as measured by IceTop. 
The plots show data at L3 level, applying the same cut 
values as in the IC40 analysis. The left plot shows the 
distribution of SLC hits for all events, while the right 
plot shows the same distribution but restricted to the 
subset of events which contain only SLC hits, i.e. events 
with no HLC hits. Hits associated with the muon bundle 
are seen throughout the detector, although the number of 
hits varies with depth because of variations in the optical 
properties of the ice due to naturally varying levels of 
contaminants such as dust, which attenuate Cherenkov 
photons The large number of isolated hits in the two 
bottom rows is an edge effect: the DOMs have fewer 
neighbors, so the chance for a hit to be isolated increases. 
In principle, the same effect could occur at the top two 
rows. However, the muon bundle deposits more energy 



in this region and the probability for any hit to have 
neighbor hits is larger here. 

The muon-poor showers that produce no HLC hits 
(right-hand plot) can still cause some isolated hits in the 
top of the detector. These events can be removed with 
an additional cut on SLC hits. Because isolated hits can 
also be produced by noise, only a small area is selected in 
which SLC hits are used as a veto. A simple additional 
L4 cut is that all events are removed that have an SLC 
hit meeting the following three criteria: 

• it is within 200 m from the reconstructed shower 
axis, 

• it is within a time window of 4.8-7.5 /is after the 
shower arrival time, and 

• it is in one of the six top layers of DOMs (spanning 
a vertical extend of 85 m). 

Note that the lower bound of the time window (4.8 /is) 
corresponds to the time it takes for a muon traveling 
vertically to reach the top layer of in-ice DOMs starting 
from the surface. Muons from an inclined shower will 
arrive even later. The number of background events that 
are discarded in the L4 cut is increased by ~30%, while 
the SLC noise rate in the data implies a decrease in signal 
efficiency of ~5%. 

With the completed detector it will be possible to op- 
timize the SLC cuts further by making the time window 
dependent on zenith angle and DOM depth. The effect 
of this optimization was not yet studied here. 



C. Re-optimization of cuts 

For the IC40 analysis the cut parameters were opti- 
mized to increase the detection probability of a possi- 
ble correlation of gamma rays with the Galactic Plane. 
To increase the sensitivity of future searches with the 
completed IC86 configuration, the cut values were re- 
evaluated to increase the number of candidate events 
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FIG. 11. Observed time of isolated hits (SLC) relative to the arrival time of the air shower front measured by IceTop. The 
left plot shows the distribution of SLC hits for all events; the right plot is the same but restricted to the subset of events which 
have only SLC hits and no HLC hits. There is an excess of SLC hits in the region were a muon signal associated to the shower 
is expected. This allows an additional cut to separate gamma-ray showers from hadronic showers. The black box indicates the 
region in which an SLC hit counts as a veto (see text for details). The variation in the number of hits as a function of depth 
in the left plot is due to variations in the optical properties of the ice. 




Energy (GeV) Distance from inner string (m) Icetop Containment 

FIG. 12. The number of L3 (red) and L4 (blue) events in the data as a function of the three main cut parameters. The ratio 
of the number of L4 to L3 events is given by the black dotted line, for which the corresponding axis is drawn on the right-hand 
side of the plot. 



without losing background rejection power. This was 
achieved by evaluating the ratio between the number of 
events after L3 and L4 cuts. 

While the L3 event sample is completely dominated 
by CRs, the L4 sample is a combination of possible 
gamma-ray showers, muon-poor CR showers and mis- 
reconstructed CR showers. The fraction of gamma-rays 
and muon-poor CRs in the detected events is independent 
from the cuts on geometry parameters d str and Crr- The 
number of misreconstructed CR showers, on the other 
hand, will increase if the geometry cut values are cho- 
sen too loosely. Therefore, the ratio between the number 
of L4 and L3 events as a function of the cut parameter 
should be fiat up to some maximum value after which it 
starts to increase. This maximum value is the preferred 
cut value since it maximizes the number of candidate 
events without lowering the background rejection power. 
It also maximizes the FOV, as looser geometry cuts imply 
a larger maximum zenith angle. 



Figure 12 shows the number of L3 (red) and L4 (blue) 



events together with their ratio (black dotted line; right- 
hand axis) as a function of the three main cut parameters 
(with the other cut parameters kept constant at their 
final value). 

The rejection efficiency for d str is fairly stable up to 
60 m. The number of events rapidly decreases above this 
value, while the rejection becomes worse. In this case, 
the alternative method of optimization yields the same 
result as the method used in the IC40 analysis. For the 
containment size Crr the ratio remains stable up to the 
edge of the array (Crr = 1) after which it starts to rise. 
It appears the cut can be relaxed with respect to the 
IC40 analysis. In the following we will use d s t Y < 60 m 
and Cit < 1.0. 

The efficiency of the energy cut increases, as expected, 
with increasing energy, leveling off around ~ 2.0 PeV. 
Since the total number of events falls off rapidly for in- 
creasing energy, the most sensitive region will be ~ 2 — 3 
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FIG. 13. Acceptance (effective area integrated over solid 
angle) for showers with an axis through both IceTop and Ice- 
Cube for IC40, Cit=0.8 (black), and IC86, Crr=1.0 (blue). 



PeV. However, since the spectra of possible sources in 
this energy regime are unknown, it is not clear what en- 
ergy cut would produce the optimal sensitivity. Instead, 
the sensitivity is calculated for ten energy bins in the 
range 1-10 PeV (see Fig. 14). 



D. Increased acceptance 

With a larger array the acceptance, defined here as 
the effective area integrated over the solid angle of each 
1° bin in zenith angle, increases considerably. Because 
of the condition that the shower axis has to be inside 
the instrumented area of both IceCube and IceTop, the 
increase is especially dramatic at larger zenith angles. 
Fig.[T3"1shows the acceptance for IC40 with Cyr < 0.8 and 
the complete IC86 array with Cn < 1.0. Not only does 
the acceptance increase at large zenith angles, the range 
of possible zenith angles is also extended (to w 45°). This 
extends the FOV to cover a larger part of the Galactic 
Plane and probe an area closer to the Galactic Center. 
The Galactic Center itself is still outside the FOV at 
5 m —29°, corresponding to a zenith angle of 61°. 



E. Sensitivity 

The sensitivity that can be reached with 5 years of data 
from the completed IceCube configuration can be esti- 
mated with preliminary data from IC79. It is assumed 
that the fraction of gamma-rays that are missed due to 
noise hits is the same as in the IC40 analysis. The full 
detector obviously has more noise hits, but this can be 
compensated by refining the in-ice cut by only allowing 
vetoes from DOMs that can be associated to the shower 
muon bundle in space and time (cf. the SLC cut described 
above). The sensitivity is calculated by producing scram- 



bled data sets with randomized RA. Figure 14 shows the 



FIG. 14. Existing limits (red triangles for CASA-MIA and 
purple line for present IC40 analysis) and IceCube sensitivity 
to a diffuse gamma-ray flux from a region within 10 degrees 
from the Galactic Plane. The blue dashed line indicates the 
five year sensitivity of the completed detector, while the blue 
dots represent the sensitivity in smaller energy bins. 



90% C.L. sensitivity to a diffuse flux from within 10 de- 
grees of the Galactic Plane that can be achieved with 5 
years of full detector data. The blue dashed line indi- 
cates the integrated limit between 1 and 10 PeV, while 
the blue dots indicate the sensitivity in six smaller en- 
ergy bins. The upper limits found by CASA-MIA and 
IC40 (present work) are also included in the plot. The 
KASCADE results are not included since they set a 
limit on the all-sky gamma-ray flux. 

Figure [l~5| shows the sensitivity to point sources that is 
possible with 5 years of IceCube data. The sensitivity is 
a strong function of declination because the acceptance 
decreases at larger zenith angles. Point sources are ex- 
pected to lie close to the Galactic Plane which reaches 
its lowest declination at —63°. Within the IceCube field 
of view there are several PWNe and other gamma-ray 
sources detected by H.E.S.S. [27], listed in Table [!] For 
these sources no significant cut-off was observed up to 
the maximum energy of 10 TeV, where statistics gets low. 
The blue dots indicate the flux that these sources would 
have at 1 PeV if their spectrum remains unchanged up 
to that energy. No correction for gamma-ray attenuation 
between the source and observer has been applied in this 
calculation. The extrapolation over two order of mag- 
nitude causes large uncertainties in the gamma-ray flux 
due to propagation of the errors on the spectral indices. 



VII. CONCLUSIONS 

We have presented a new method of searching for high 
energy gamma-rays using the IceCube detector and its 
surface array IceTop. One year of data from IC40 was 
used to perform a search for point sources and a Galac- 
tic diffuse signal. No sources were found, resulting in 
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TABLE I. List of H.E.S.S. sources in IceCube gamma-ray FOV. For those values that have two error margins, the first indicates 
the statistical error, while the second indicates the systematic error. 
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FIG. 15. IceCube 5 year sensitivity to point sources as a 
function of declination. The solid (dashed) black line indicates 
the sensitivity to an E~ 2 (E~ 2 ' 5 ) flux. The dashed red line 
indicates the lowest declination reached by the Galactic Plane. 
The blue points indicate the flux at 1 PeV with extrapolated 
uncertainties of the sources listed in Table [I] in the absence of 
a cut-off. 



a 90% C.L. upper limit on the ratio of gamma rays 
to cosmic rays of 1.2 x 1CP 3 for the flux coming from 
the Galactic Plane region ( -80° < I < -30°; -10° < 
b < 5°) in the energy range 1.2 - 6.0 PeV. Point source 
fluxes with E~ 2 spectra have been excluded at a level of 



(£/TeV) 2 d$/d£; ~ 10 
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cm-^s^TeV- 1 de- 



pending on source declination. The full detector was 
shown to be much more sensitive, because of its larger 
size, improved reconstruction techniques and the possi- 
bility to record isolated hits. 

This analysis offers interesting observation possibili- 
ties. IceCube can search for a diffuse Galactic gamma-ray 
flux with a sensitivity comparable to CASA-MIA, but at 
higher energies. This sensitivity is reached, however, by 
studying a much smaller part of the Galactic Plane than 
CASA-MIA. IceCube is therefore especially sensitive to 



localized sources, which might be Galactic accelerators 
or dense targets for extragalactic CRs. 

The H.E.S.S. and CANGAROO-III [38] telescopes 
have found several high energy gamma-ray sources in 
IceCube's FOV. Most of these sources are identified as 
or correlated with PWNe. Their energy spectrum has 
been measured up to a couple of tens of TeV. At this 
energy, statistics become low and for most sources no 
cut-off has been established. If these spectra extend to 
PeV energies without a break, IceCube will be able to 
detect them. It is also possible that an additional spec- 
tral component in the PeV energy range is present if a 
nearby dense molecular cloud acts as a target for the 
PWN beam [39]. IceCube will be able to study these 
systems and place constraints on their behavior at very 
high energies, or possibly detect PeV gamma-rays for the 
first time. 
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